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ABSTRACT 


The multiple mesh/ single stage dynamics program is a gear tooth 
analysis program which determines detailed geometry, dynamic 
loads , stresses and surface damage factors . The program can 
analyze a variety of both epicyclic and single mesh systems with 
spur or helical gear teeth including internal , external , and 
buttress tooth forms. 

The modifications refine the options for the flexible carrier and 
flexible ring gear rim and adds three new options: a floating sun 
gear option, a natural frequency option, and a finite element 
compliance formulation for helical gear teeth. The option for a 
floating sun incorporates two additional degrees of freedom at the 
sun center. The natural frequency option evaluates the frequencies 
of planetary, star or differential systems as well as the effect 
of additional springs at the sun center and those due to a 
flexible carrier and/or ring gear rim. The helical tooth pair 
finite element calculated compliance is obtained from an automated 
element breakup of the helical teeth and then is used with the 
basic gear dynamic solution and stress postprocessing routines. 
The flexible carrier or ring gear rim option for planetary and 
star spur gear systems allows the output torque per carrier and 
ring gear rim segment to vary based on the dynamic response of the 
entire system, while the total output torque remains constant. 
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ABSTRACT 


The multiple mesh/single stage dynamics program is a gear tooth 
analysis program which determines detailed geometry , dynamic 
loads, stresses and surface damage factors. The program can 
analyze a variety of both epicyclic and single mesh systems with 
spur or helical gear teeth including internal, external, and 
buttress tooth forms. 

The modifications refine the options for the flexible carrier and 
flexible ring gear rim and adds three new options: a floating sun 
gear option, a natural frequency option, and a finite element 
compliance formulation for helical gear teeth. The option for a 
floating sun incorporates two additional degrees of freedom at the 
sun center. The natural frequency option evaluates the frequencies 
of planetary, star or differential systems as well as the effect 
of additional springs at the sun center and those due to a 
flexible carrier and/or ring gear rim. The helical tooth pair 
finite element calculated compliance is obtained from an automated 
element breakup of the helical teeth and then is used with the 
basic gear dynamic solution and stress postprocessing routines. 
The flexible carrier or ring gear rim option for planetary and 
star spur gear systems allows the output torque per carrier and 
ring gear rim segment to vary based on the dynamic response of the 
entire system, while the total output torque remains constant. 
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I . SUMMARY 


The epicycllc gear program is a multiple mesh/single stage, gear 
dynamics program. It is a versatile gear tooth dynamic analysis 
computer program which determines detailed geometry, dynamic 
loads, stresses and surface damage factors. The program can 
analyze a variety of both epicyclic and single mesh systems with 
spur and helical gear teeth including internal, external, and 
buttress tooth forms. This NASA Lewis sponsored contract called 
for four improvements: refinement of the option for the flexible 
carrier or flexible ring gear rim, a floating sun gear option, a 
natural frequencies option, and a finite element compliance 
formulation for helical gear teeth. 

Task I was to add an option for a floating sun to account for 
flexible sun gear mounting , which incorporates two additional 
degrees of freedom at the sun center. Generally, soft mounted sun 
gears are used to minimize the effects of gear runout, etc. The 
test case used for the program checkout was a lightly loaded, 
three planet, planetary system. The floating sun case results 
showed similar loads for the sun-planet meshes and slightly higher 
loads for the ring-planet meshes when compared to the rigidly 
mounted sun case . Other cases , with higher loads and various 
spring rates, should be examined to analyze the effect more 
thoroughly . 

Task II was to add an option to determine the natural frequencies 
of the system. This option is desirable to predict system 
critical speeds without having to run the complete dynamic 
analysis for a range of input speeds. The information can also aid 
the user in running the dynamic solution routines. At the system 
critical speeds the dynamic loads will be much more sensitive to 
the input variables and the user will know where these speeds 
occur beforehand. Planetary, star or differential systems can be 
investigated as well as the effect of the additional springs at 
the sun center (floating sun) and those due to a flexible carrier 
and/or ring gear rim. In addition, the effect of variation of the 
tooth pair stiffness on the natural frequencies due to load 
position can be investigated. The frequency results are 
consistent with the program's dynamic response solution. 

Task III was to generate a helical tooth compliance routine based 
on finite element modeling. The previous version divided the 
tooth into ten equivalent spur gear teeth and used the spur tooth 
routines for the dynamic solution. However, this technique did 
not allow for coupling between the equivalent spur teeth segments. 
The program also had provisions to input a general compliance 
matrix for the helical gear tooth to be analyzed, but the user had 
to know the matrix before running the gear dynamics program. This 
added finite element routine eliminates these prior shortcomings. 
The new option internally generates a finite element breakup of 
the helical teeth and the necessary data for the internal finite 
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element routines. The routines use a four noded, quadrilateral, 
higher order plate element with five degrees of freedom per node. 
The results are used to obtain a general tooth pair compliance 
curve which is then used by the basic dynamic solution routines . 

Task IV was to refine the flexible carrier and ring gear rim 
options for planetary and star spur gear systems. The frequency 
equations were expanded to account for nonrigid carriers and ring 
gear rims. In addition, some minor modifications were made with 
respect to the numerical solution tolerances. These modifications 
results in more stable solutions, thus allowing the user to 
investigate the effects of various spring rates for both the 
carrier and the carrier /planet pin, eg. a bearing, or for the ring 
gear rim and the corresponding pin. 


II. 


INTRODUCTION 


A. PROGRAM HISTORY 

The multiple mesh gear dynamic analysis computer code has been 
under development at Hamilton Standard for about five years. The 
program can determine detailed geometry, dynamic loads, stresses, 
and surface damage factors for epicyclic gear systems and single 
mesh systems with internal, external, buttress, or helical tooth 
forms. The significant parameters can be plotted through the 
entire mesh in addition to the maximum values which are tabulated 
as output from the program. 

The initial program, a single spur gear mesh, was written for high 
contact as well as low contact ratio gearing. The basic concept 
was an extension of ' that developed by Richardson in 1958 . Since 
the basic program was developed, many enhancements and refinements 
have been made. Buttress, internal and external involute tooth 
forms can be analyzed for spur and helical gear teeth. Tooth 
spacing errors, runout errors, involute profile modifications, 
etc. can be accounted for in the dynamic solution. 

The program has been developed to operate over a wide range of 
contact ratios, and to allow the gear teeth to have a different 
pressure angle on the drive side and coast side (buttress form) . 
The teeth of the meshes being analyzed may have modified profiles 
and spacing errors may be specified. Influence coefficients may be 
input for rim and web geometry to determine the effect on peak 
dynamic stress caused by non-uniform stiffness along the width of 
the gear, dependent on the design configuration of the foundation 
under the gear tooth. More recent additions to the program include 
variable contact friction throughout each mesh, user friendly 
options, dynamic side bands, a speed survey option and the option 
of solving non-planetary or single mesh systems. See References 1 
through 6 for more details. 


B . PROGRAM ENHANCEMENTS 

This NASA contract refines the option for the flexible carrier or 
flexible ring gear rim and adds three new options: a floating sun 
gear option, a natural frequency option, and a finite element 
compliance formulation for helical gear teeth. 

Task I was to add an option to allow the sun to float at the 
center. This allows for investigation of the effects of different 
spring rates and damping at the sun center on the dynamic tooth 
load behavior. This was accomplished by adding two global 
translational degrees of freedom at the sun center . These 
translational degrees of freedom were transformed into degrees of 
freedom along the respective meshing lines of action, to remain 
consistent with the existing code. The theoretical development was 
aided by the work of Hidaka et al. Reference 7. 
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Task II was to develop a routine for calculating the natural 
frequencies of the gear system. This was accomplished by solving 
the classical eigenvalue problem. This option will solve for the 
frequencies of the system using tooth pair compliances at the 
pitch diameter. The designer can also utilize the load 
position varying gear tooth compliance formulation and evaluate 
the frequencies at other mesh positions. This option allows the 
user to investigate the effects of different spring rates and 
masses throughout the gear system on the critical speeds. 

Task III was to refine the helical gear tooth compliance routines. 
The refinement was to incorporate convective (coupled) compliance 
effects which were not previously included. This involved adding 
routines to build finite element models of the helical gear teeth, 
Reference 8. The models are used to obtain a spur gear compliance 
formulation which utilizes the basic gear dynamic routines and 
stress postprocessors. This refined approach uses much less CPU 
time, as well as using a smaller time step for the numerical 
solution, which will lead to a more stable dynamic solution than 
the previous uncoupled spur tooth segment approach. However, for 
helical gears with large helix angles the stress postprocessing 
will give unconservative results . 

Task IV was to refine the flexible carrier /ring gear rim option. 
This option allows the user to investigate the effects of various 
stiffnesses for the carrier or ring gear rim as well as the pin 
stiffness between the planet gear and the carrier or the ring gear 
rim and the output shaft. This involved investigation of previous 
work. Reference 1, in order to determine the cause of an 
instability. The numerical solution technique was evaluated, as 
well as review of the mathematical model. The most significant 
change was in the mathematical model, where the forcing functions 
of the carrier or ring gear segment equations were modified to 
vary with respect to time while the total output torque remained 
constant. In addition, some minor changes were made to the 
numerical solution parameters to increase stability. 

These four improvements further enhance the flexibility of the 
multiple mesh gear dynamics program and the variety of 
applications that can be modeled. The modifications allow the 
user to investigate a wider variety of system complexities. For 
example, a simple planetary can now be evaluated utilizing various 
additional degrees of freedom such as the floating sun gear or the 
flexible carrier, which are likely to be influential in a real 
system. The finite element compliance formulation lays the 
groundwork for a more exact modeling of the helical gear teeth. 
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III. FLOATING SUN GEAR 


A. PROGRAM MODIFICATIONS 

The floating sun gear option adds two additional degrees of 
freedom at the sun center. Springrates are required for the two 
orthogonal directions at the sun center as well as the 
translational mass of the sun gear and the additional boundary 
conditions. The equations of motion for the sun center are 
included in the system of equations being solved. 

The equations of motion for the sun gear center in Cartesian 
coodinates, x and y directions, using the model of Figure 1, can 
be written: 


m x + d x + k x - £ d sin a, 

s x x i . 1 sp i i sp i i 


Cl) 


- I L sin a, 
i-1 s *t 1 


“ sY + V + V + ^ d s Pl *s Pi =°s 


+ E L ST3 C0S “4-0 

i-1 S Pi 1 * 


( 2 ) 


where 


a i ■ ®o + ' r i - * 

a i - *1 - * - e r 


(planetary) 

(star) 


The angle, a., is for resolution of the tooth pair forces from 
along the line of action to the Cartesian coordinates, x and y. 
The pressure angle, §, was assumed to remain constant through the 
mesh; however, the more precise formulation would account for the 
varying angle at different mesh times. It is believed this 
assumption has secondary effects on the results. However, the 
information is available in the code such that the varying angle 
could be included during future enhancements. 

In addition to the sun center equations , the carrier or ring 
rotational displacements, corresponding to a planetary or star 
system respectively, are also required to obtain the angle, a. 
Thus the following equations must also be solved for a rigid 
carrier or rigid ring gear rim in conjunction with the tooth pair 
mesh equations. Reference 2, and equations (1) and (2). 
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( 3 ) 


(4) 


Y r - »b r ®P 

The numerical solution requires these equations be reduced to 
first order differential equations; thus equations (1), (2), and 
(3) or (4) are added to the system of equations via six first 
order equations, see Appendix B. 

The sun center displacements must be resolved in the direction of 
the line of action to determine the effect on the tooth pair 
meshes. This is accomplished via: 


X LOA t ’ J 008 “i " x Sin a i 

The tooth pair meshing loads of Reference 2, then become: 


Lsp l ' 3 -l t(Ysp i ' ^Pji " * sPji 6 sPji + 


Lrp i ' >x [(Yrp i " * rp ji ‘ * a *Pji 


(6a) 

(6b) 


where m is the number of teeth in contact at mesh i, ti is the 
tooth pair stiffness and is a tooth pair contact identity 
function . *i 

To obtain a steady state solution, the solution is iterated until 
the boundary conditions converge for the tooth pair meshes . 
Convergence is determined by comparing the displacements along the 
lines of action due to the sun center movement to the largest 
sun-planet tooth pair displacement. Convergence is faster when the 
spring rates at the sun center are of similar order of magnitude 
to the tooth pair stiffnesses. 
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B. DISCUSSION OF FLOATING SUN RESULTS 


The addition of the floating sun gear option allows the user to 
analyze various sun center spring rates and damping and the 
resulting dynamic tooth loads. Most epicyclic gear systems have a 
sun gear that can move in the in-plane translational directions, 
thus this option incorporates degrees of freedom that are of 
practical interest. 

Two test cases were run using a lightly loaded planetary gear 
system. The description of the spur gear test case, Task I, 
Example 1.1, is described in Table 1 and was run using both the 
three planets of the system and reducing the number of planets to 
two. Several spring rates at the sun center were examined. The 
results for the three planet cases are summarized in Table 2. 


The three planet cases included phasing constants to account for 
the different location of each planet mesh on its respective line 
of action. Because of the interactive dynamics, the total 
distances moved along the lines of action, due to the additional 
sun center movement, can be different for each mesh. In the three 
planet case with sun center springs 4.4 times the sun-planet tooth 
pair stiffness (10,000,000 lb. /in. and damping of 5 %) the maximum 
loads increased from 2 to 7 percent for the sun-planet loads and 
from 14 to 15% for the ring-planet meshes. The ring-planet meshes 
also indicated a decrease in the dynamic contact ratio, which is 
calculated by determining the total tooth contact time for one 
tooth pass during the dynamic solution. A low dynamic contact 
ratio represents tooth pair separation, which indicates tooth 
bouncing may be occurring, possibly causing the load increase. 
The additional degrees of freedom at the sun center may have 
introduced natural frequencies near the operating speed which 
could also cause the load increase. 

The sun center spring rates were increased to 30,000,000 lb/ in, or 
13 times the sun-planet tooth pair stiffness. The no tooth error 
case showed the sun-planet loads less than 3 % different and the 
ring-planet loads were generally a little higher, 4-5 %, for 2 % 
damping. Thus, as the springs become stiffer, the loads approach 
the non-flexible sun center mount solution as expected. 

The two planet case showed no variation in the maximum tooth pair 
loads. As expected, the diametrally opposed planets (equal 
phasing constants) moved in equal and opposite directions along 
their lines of action due to equal tooth pair stiffnesses. 
The two cases executed to confirm these results were for sun 
center spring rates of 300,000 lb/in and 30,000 lb/in. Both cases 
yielded results identical to the non-floating sun gear. 
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IV. NATURAL FREQUENCY OPTION 


A . PROGRAM MODIFICATIONS 

The natural frequency option allows the user to investigate 
natural frequencies of the system through a classical eigenvalue 
solution. The effects of various spring rates at the sun center 
for a floating sun, or the effect of various stiffnesses of a 
planet carrier and/or ring gear rim on the frequency response can 
readily be determined. In addition, by using the planet phasing 
constants , the range of frequencies due to the nonlinear tooth 
meshing action can be investigated. 

The general form of the dynamic equations for the eigenvalue 
solution is: 


t m + mrt 


(7) 


The mass and stiffness matrices [M] and [K] , are derived from the 
system of equations of References 1 and 2, as well as equations 
(1) to (4). The eigenvector /eigenvalue solution, which assumes 
harmonic motion, then solves for the roots of the determinant: 


[ K ] - X [ M ] 


( 8 ) 


A standard eigenvalue/eigenvector numerical solution routine for 
real, symmetric matrices was used to solve the determinant. 
Reference 9. 

The program also calculates the gear mesh frequencies using the 
following equation. 


fj - rpm/60. * XN • I X-l, 2, 3... 


These gear mesh frequencies are printed with the natural 
frequencies to aid the user in generating critical speed diagrams. 

The natural frequency option can be used for eight system types: 
planetary, star, differential, single meshes, planetary with 
flexible carrier, star with flexible ring gear rim, or a 
differential system with both a flexible ring gear rim and a 
flexible planet carrier. In addition, the floating sun degrees of 
freedom can be included in the natural frequency solution. 


This option is initiated by a trigger with two choices of output 
format, either natural frequencies and gear mesh frequencies or 
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full output which also includes the eigenvectors. After the 
frequencies have been calculated, the program ends and does not 
continue with the dynamic load solution. 

The frequency solution uses tooth pair stiffnesses from the 
nonlinear compliance formulation of the gear program. Reference 4 
and 5 . This compliance formulation also includes the Hertzian 
effect; thus, the torque that is input will influence the natural 
frequencies. A zero torque case will eliminate the Hertzian 
effect if it is not desired. 

The natural frequency solution is for a specific instant in time, 
with the corresponding tooth pair stiffnesses. If no phasing 
constants are included, the tooth pair stiffnesses will be those 
at the pitch diameter. The phasing constants can be used to 
simulate different times and therefore different stiffnesses in 
the mesh. These phasing constants tell the program that the 
different planet meshes are at different positions with respect to 
each other along their lines of action. This indicates that all 
the planet meshes may not be at the pitch diameter initially. See 
User's Manual for details on calculation of the phasing constants. 


B. DISCUSSION OF NATURAL FREQUENCY RESULTS 

The natural frequency option is a useful tool for predicting 
critical speeds . This option allows for investigation of the 
effect of mesh position dependent tooth pair stiffnesses on the 
natural frequencies, as well as the additional frequencies due to 
carrier and ring flexibilities and floating sun flexibilities. 
The eigenvalue /eigenvector solution executes rapidly, as the 
program does not continue with the dynamic load solution. Thus, 
this is an economical approach to investigate the effects of 
various spring rates or masses on the natural frequencies. 

The natural frequency option can be used to calculate the speed 
ranges where high dynamic loads could occur. This can also assist 
the user in reducing the number of boundary condition iterations 
that are necessary for convergence for a regular dynamic load 
solution by avoiding these critical speed areas. The effects of 
different spring rates throughout the gear system on the critical 
speeds can also be investigated. It also allows the designer to 
increase the number of degrees of freedom and observe the 
additional frequencies. For example, a simple planetary system 
will have N + 2 degrees of freedom, where N is the number of 
planets . The designer can add two spring rates at the sun 
center for N + 4 degrees of freedom, then could increase the total 
degrees of freedom again to 2N + 4 via the carrier flexibility. 

Several test cases were run, and two examples are presented 
in Figures 3 and 4. Both figures show the critical speeds 
predicted by the natural frequency option, and the critical speeds 
predicted by running speed surveys with the dynamic response 
solution. Reference 1. Both cases agree in predicting the natural 
frequencies when the pitch diameter tooth pair stiffnesses were 
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used for the eigenvalue solution. 


Figure 3 shows a critical speed diagram for Example 2.1, Table 1, 
and Figure 4 shows Example 2.2. The horizontal bands illustrate 
the range of frequencies that result from the variation of tooth 
pair stiffnesses through the mesh. The variation in frequency due 
to tooth pair stiffnesses versus percent of total mesh positions 
is shown in Figure 5 and corresponds to Figure 3. The higher 
frequencies of the band are when the mesh is at the pitch diameter 
and the tooth pair stiffness is maximum. The lower frequencies of 
the bands are at other positions in the mesh. These were simulated 
by adjusting the planet phasing constants until the minimum 
stiffness for the sun-planet meshes and the minimum stiffness for 
the ring-planet meshes were obtained. 

The case in Figure 3 showed reasonable correlation between the 
natural frequency option and the results of an analytical speed 
survey using the dynamic solution. This was a two planet case 
with unequal phasing. The speed survey response results indicated 
peat dynamic loads in the range of sun gear input speeds of 23 , 000 
to 25,000 rpm, a distinct peak near 7,000 rpm and a smaller load 
increase near 12,000 rpm. This case was previously shown to have 
reasonable correlation to test data, Reference 1. 

By varying the stiffnesses the frequencies can be easily 
associated with particular degrees of freedom; e.g., in Example 
2.2, when the sun-planet stiffness decreased significantly, the 
highest frequency decreased and when the ring-planet stiffness 
decreased, the two lower frequencies decreased. Thus, the highest 
of the gear mesh frequencies is due to the sun-planet meshes , 
while the lower frequencies are due to the ring-planet meshes . 
This may not be immediately obvious from the eigenvectors, because 
they indicate the motion along the lines of action for each gear 
and not the relative motion of the sun-planet or ring-planet 
meshes . 

For planetary, star, and differential systems, the carrier and/or 
ring gear are treated as rigid bodies, unless the flexible carrier 
or flexible ring gear rim options are selected. The rigidity is 
also evident in the resulting natural frequencies. The eigenvalue 
solution yields one rigid body mode for planetary and star 
systems corresponding to the rigid carrier or ring gear, and the 
differential system results in two rigid body modes corresponding 
to both the rigid carrier and ring gear. These rigid body modes 
are not shown in the table of natural frequencies; however, if 
the user requests full output, all of the frequencies are printed 
as well as the mass and stiffness matrices and all eigenvectors. 


Systems with a flexible carrier and/or ring gear rim may yield 
rigid body modes if the stiffnesses of the carrier and/or ring 
gear rim or pin stiffnesses are relatively high. For these 
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situations, either the carrier and/or the ring are acting as rigid 
bodies. For these systems all the frequencies are output, but if 
the first two modes are orders of magnitude less than the other 
frequencies, the pin stiffnesses and/or the program results should 
be carefully examined. 


V. REFINED HELICAL COMPLIANCE ROUTINES 


A . PROGRAM MODIFICATIONS 

An option for helical gear analysis was added which incorporates a 
finite element analysis to obtain a tooth pair compliance curve. 
Several new subroutines have been added to the code to generate 
finite element models of the helical gear teeth. These models are 
then used to calculate tooth pair compliance curves. The model 
generation is internal to the code, so the user need only input an 
additional trigger to initiate the option. 

The finite element option uses two finite element plate routines, 
Reference 8, to generate the stiffnesses for both in-plane and 
out-of-plane loads, as well as a routine to process the 
information to and from the multiple mesh program . Figure 6 
summarizes the procedure via a flow chart. 

For both routines, four noded isoparametric plate elements, with a 
total of five degrees of freedom per node, are used. The routine 
for out-of -plane loads includes transverse shear effects. The 
model for any tooth will have 9 elements along the face width and 
9 elements along the tooth centerline (10 by 10 nodes), see Figure 
7. The plate thicknesses are average thicknesses determined from 
the existing tooth geometry subroutines. 

The tooth model is fixed at the root and displacements are applied 
along seven equally spaced load lines via boundary conditions for 
the finite element solution. The reaction forces along the load 
lines are used in conjunction with the applied displacements to 
obtain average stiffnesses for seven load line positions. 

The existing spur tooth pair compliance used in the program 
includes axial bending, Hertzian, and fillet and foundation 
effects. Reference S and 6. These are determined for seven load 
positions for each tooth pair and combined to obtain total 
compliances. A curve fitting routine is then used to obtain the 
following fourth order polynomial for compliance as a function of 
position along the line of action. 


C - C 0 [ 1 + A(S/S 0 ) + B(S/S 0 ) 2 + C(S/S 0 ) 3 + D(S/S 0 > 4 1 (10) 


The finite element model accounts for the axial bending and 
fillet effects and an average Hertzian effect. The load line 
position at the center of the helical tooth face width is used to 
calculate an average Hertzian compliance. 
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This approach accounts for the helix angle, but uses the spur gear 
dynamic solution technique. The previous helical solution divided 
the tooth into ten equivalent spur gear tooth segments. Each of 
the segments was evaluated dynamically and the CPU running time 
was quite long to solve for each segment at each of 100 time 
steps. Because of this, the number of time steps had been reduced 
to ten for the previous helical solution. However, because the 
finite element approach is more direct than the segment approach, 
the time step was increased back to 100 for this option to improve 
the accuracey and response definition. 


B. DISCUSSION OF REFINED HELICAL GEAR COMPLIANCE RESULTS 


This option offers the user an alternative helical gear tooth 
analysis. A compliance curve was generated using the existing 
formulations for Hertzian effects and detailed geometry to obtain 
plate thicknesses to be used with two finite element routines , one 
for in-plane loads and one for out-of-plane loads, see Figure 8 
for example output of Task III, Example 3.1, Table 1. 

Future improvements should include further refinement of the 
stress postprocessing, where the finite element routines could be 
utilized more fully. The current finite element configuration 
utilizes the spur gear stress postprocessing routines by using the 
calculated dynamic load and applying it to the center of the tooth 
face width. This could lead to unconservative stress results when 
the helix angle is large, because the load line is parallel to the 
axis of rotation for spur gear postprocessing. A more complete 
solution would utilize the flexibility matrix generated by the 
finite element routines in conjunction with the dynamic solution. 
This more precise method would involve calling the finite element 
routines directly during the dynamic solution, which would require 
significant amounts of additional computer time. Similarly, 
finite element stress sensitivity routines that use the plate 
element directly for postprocessing could be included. 

The number of elements that are chosen for the mesh breakup should 
be changed to be a user input in future work. The current constant 
breakup could lead to width to height ratios that exceed 
acceptable limits for finite element aspect ratios for wide teeth. 
The tooth fillet elements could also be refined. It was assumed 
that the first row of elements would extend from the root to the 
first point of contact, while the plate thicknesses were modified 
to estimate the additional fillet material, see Figure 9. However, 
this does not accurately define this area and could lead to 
significant element size variation and therefore error for cases 
with large fillet radii relative to the tooth height. 
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VI. FLEXIBLE CARRIER EVALUATION 


A. PROGRAM MODIFICATIONS 

During the previous contract work (CR # 174747) an attempt had 
been made to add the flexible carrier option to the dynamics 
model. However, numerical instabilities had occurred in the check 
case solutions which were not resolved during that contract . 
Therefore, the flexible carrier/ring gear rim dynamic solutions 
were reviewed with respect to the numerical solution, the 
equations, and the FORTRAN code. The major improvement was to the 
output torque term of the flexible carrier or flexible ring gear 
equations, i.e. the forcing function for the differential 
equation. This term was modified to calculate the output torque 
from the carrier or ring segment displacements and an interface 
stiffness, such as a pin or bearing stiffness, while the total 
output torque was constrained to remain constant. 

The theoretical model is shown in Figure 10, and can be used to 
write the carrier or ring segment equations of motion as follows. 


VV W ’■ ' V V - 'VV v v 

- l «,- s t * vv v,’ 
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Equation (13) adds one more equation than there are unknowns; 
thus, a constraint must be simulated in the numerical solution 
routines by introducing another parameter. The numerical routine 
used to solve the system of first order equations does not have 
the direct capability of a constraint on a primary variable, i.e. 
Y C 4* Therefore, an artificial mass which simulates a large 
inertia, assumed to be 1000 times the largest mass of the system 
and attached to the carrier ^ or ring, was added to the system of 
equations. This results in x - 0, analogous to the system shown 
in Figure 12. The additional equation of motion is: 
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N 
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or: 


(15) 

(ring) 
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This increases the number of second order equations by 1 to (3N + 
2). The first order equations used for the numerical solution 
are in Appendix B. 

The numerical solution routine was also investigated with respect 
to stability. A simple example was used in conjunction with the 
numerical solution routines to determine the effect of the size of 
the integration time steps on the solution. The example used a 
system of differential equations similar to the equations for a 
planetary system. They were used directly in the numerical 
solution routines, independent of the gear program. It was 
verified that the time step used for the numerical integration 
could lead to erroneous results if the step size was too large. 
This step size was previously a constant value in the multiple 
mesh code, but the total mesh time varies for different systems. 
A better choice for the integration time step is a function of the 
meshing time. Therefore, the integration time step was then 
changed to be 0.001 percent of the time for each mesh cycle. 

This step size change can actually reduce the number of iterations 
internal to the numerical routines, due to the way the solution 
routine handles the step size if it is too large. If the routine 
is given a step size that is too large, it reduces the integration 
step size and tries again until it is small enough to yield a good 
solution. If it is given a step size that is adequate initially, 
it will obtain a good solution more rapidly. As this occurs for 
each of the 100 time steps in the dynamic solution, this could 
lead to significant reductions in computation time for some cases. 

The size of the mesh time step can also have an effect on the 
solution. This is currently set to be equal to the total mesh 
time divided by 100. For very low speeds this time step can become 
relatively large and could lead to an unstable solution. 


B. DISCUSSION OF FLEXIBLE CARRIER RESULTS 

The flexible carrier modifications allow the user to investigate 
the effect of a carrier with additional flexibility, along the 
line of action, on the dynamic tooth loads of a spur gear system. 
Similarly, additional flexiblity in the ring gear rim along the 
line of action. Two test cases with different carrier stiffnesses 
for a planetary spur system (flexible carrier) were investigated. 
They indicated higher maximum loads for the more flexible systems. 
This is due to the increased tooth pair bouncing that occurred, as 
was indicated by lower dynamic contact ratios. 
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Figures 12a and 12b illustrate the non-flexible test case 
numerical results and a case , Example 4.1, Task IV , Table 1 , with 
a pin stiffness of 2,500,000 lbs/in. and a carrier stiffness of 
5,000,000 lbs/in. The maximum tooth loads increased with the 
increased flexibility of the planet carrier. For this case, the 
ring-planet loads increased more than the sun-planet loads, 12.5% 
and 5.9%. Results from another case, with equal carrier and pin 
spring rates showed that the sun-planet mesh loads increased 9.3%, 
while the ring-planet loads increased by 5.4%. 

Figures 13a and 13b show the plots generated by the multiple mesh 
program for the first planet mesh for the nonflexible case and for 
the case with pin stiffness of 2,500,000 lbs/in., respectively. 
There were no planet phasing constants for this test case; 
therefore, all three planet mesh results are Identical and only 
the results for one planet are shown. Comparison of the flexible 
and nonflexible load plots indicated that while the maximum loads 
increased, the minimum loads remained nearly the same. There was a 
similar effeot observed in the other plotted parameters; that is, 
the curves maintained the same general shape, the minimum points 
did not noticably change, but the maximums increased with the 
addition of the flexible carrier. 

Another test case, Example 4.2, Table 1, was executed for a 
lightly loaded planetary system with unequal phasing and with a 
carrier stiffness of 3,000,000 lbs/in and a pin stiffness of 
5,000,000 lbs/in. The results for this case, Figure 12c, show 
very little variation from the rigid carrier results, less than 
2 %. This variation is most likely due to the very small loads. 

It should be noted that this solution consumes considerably more 
CPU time than a simple planetary or star system due to the larger 
system of dynamic equations and particularly the increased number 
of boundary conditions that must converge. It is recommended that 
the nonflexible solution be run first and the converged boundary 
conditions be used in conjunction with the output torque term to 
estimate the new boundary conditions, see Appendix A for details. 

Due to funding limitations, a flexible ring rim case was not 
executed to completion. It has been coded and a test of the first 
few boundary condition iterations indicated that there should be 
no significant problems with this code. 
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VII . CONCLUDING REMARKS 


A number of new options have been added to the multiple mesh gear 
program, but only a few test cases were run for each of these due 
to the minimum amount of epicyclio data available for comparision. 
It is recommended that the program be more thoroughly tested and 
evaluated via parametric studies as well as comparison to test 
data. Follow on work should also include refinement of the finite 
element stress sensitivity postprocessing. 

A floating sun gear option was added which will allow the user to 
investigate the effects of various spring rates and damping at the 
sun center . Several test cases were executed and the program 
gives reasonable results, with the solution approaching the rigid 
mount solution as the spring rates become 10 to 15 times the tooth 
pair stiffnesses. 

The critical speeds predicted by the new natural frequency option, 
and the critical speeds indicted by running speed surveys with the 
dynamic response solution agreed well with the frequencies for 
tooth pair stiffnesses at the pitch diameter. 

The refined helical gear option lays the foundation for an 
alternative method of analyzing helical gear dynamics. Due to 
funding limitations, the potential for the finite element solution 
was not fully accomplished. However, the finite element results 
were used to generate a general tooth pair compliance similar to 
the original spur gear compliance formulation. The stress 
postprocessing uses the spur gear stress postprocesser which 
inherently assumes the load line is parallel to the axis of 
rotation. Thus, the stress postprocessing could have significant 
errors for large helix angles. 

The flexible carrier and flexible ring gear rim options modified 
the output torque for each segment to vary with the dynamic 
results while the total output torque remains constant. The test 
case results were reasonable, with the maximum loads increasing 
and the minimum loads remaining nearly the same when compared to 
the rigid carrier results. 
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Figure 2. Flowchart for Natural Frequency Modifications 


Subroutine FREQ is called from READ2 
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, Figure 3. Natural Frequency Interference Plots, Example 1 
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Natural Frequency (KHz) 


Figure 4. Natural Frequency Interference Plots, Example 2 
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Figure 6. Flowchart for Helical Tooth Compliance Modifications 
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Abbreviations 


FE = Finite Element 
IP = In-Plane 
OOP = Out-of-Plane 
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Figure 7. Helical Gear Tooth Finite Element Model 
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Figure 8 . Program Output from Finite Element 

Helical Gear Option 
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Figure 9 . Fillet Element Thickness 


45° for fillet element 
thickness calculations 



R r - root radius 

Rg = base radius 

Rp = pitch radius 

R^ = tip radius 

r = fillet radius 

L = load 

0 - pressure angle 

0 # = load line angularity 

t = fillet element thickness 













Figure 11. Torque Constraint Model Analogy 
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MAXIHUM CARRIER DISPLACEMENT 0.A254E-03 


Figure 12a (cont) 
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Figure 13a. Non-flexible Carrier Gear Mesh Plots 
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0.0 distance along line of action 



sun-planet mesh 


sun-planet mesh 




Figure 13a (cont) 
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Figure 13b (cont) 
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APPENDIX A: USER'S MANUAL 


Part 1 : Summary of New Option Information 


The new options can he easily added to any existing data set. The 
new inputs corresponding to the new options are summarized below. 
An updated version of the User's Manual of Reference 1 follows in 
Part 2 with all the input descriptions. In addition, some 
comments are included that may be useful for input or 
interpretation of output. It should be noted that the locations to 
request plots have been changed. 


Floating Sun: 


Springrates (locations 115 and 116) and damping (locations 117 and 
118) are input for two Cartesian directions, horizontal and 
vertical or x and y , at the sun center . The translational mass 
of the sun gear (location 119) and the additional boundary 
conditions (locations 561 to 644) are also inputs. The additional 
boundary conditions that must be input are for the sun center and 
the carrier or ring displacements and velocities, see the updated 
location listing , Part 2 , for details on input . 


Note that some damping should be input or numerical instabilities 
may develop. 


The convergence criteria should not initially be too stringent , 
especially if the spring rates at the sun center are soft. This 
can lead to diverging boundary conditions, especially for systems 
with unequal phasing constants. 


Natural Frequencies: 

The user can utilize existing input data sets by simply adding a 
trigger (location 805). In general, minimum input (Level I) is 
all that is necessary; however, any existing data set can be used 
with the addition of the trigger. The trigger is also used to 
indicate the type of output -- frequencies only or 
eigenvalues/eigenvectors , etc. After the frequencies have been 
calculated, the program ends and does not continue with the 
dynamic solution. 


A list of the various system types that can be investigated and 
the corresponding number of degrees of freedom follows . 
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k - type of epicyolic spur gear system 

n = number of planets 

ndf * number of degrees of freedom 


k - 1 planetary system, i.e. ring gear fixed and rigid 
planet carrier 
ndf - n + 2 

k = 2 star system, i.e. planet carrier fixed and rigid ring 
gear rim 
ndf = n + 2 

k = 3 differential system, i.e. rigid ring gear rim and 
rigid planet carrier 
ndf - n + 3 

k = 4 single mesh, sun-planet (s) or external /external 
mesh 

ndf = n + 1 

k = 5 single mesh, ring-planet (s) or external /internal 
ndf - n + 1 

k = 6 planetary system with flexible carrier 
ndf - 2*n + 1 

k - 7 star system with flexible ring gear 
ndf = 2*n + 1 

k = 8 * differential system with both flexibar and 

flexible planet carrier 
ndf - 3*n + 1 

♦NOTE: For k-8, a differential system with both a flexible 

carrier and a flexible ring, there is only a natural frequency 
solution, i.e. there is no dynamic solution. 

To include the floating sun option or for systems with k of 6, 7, 
or 8, note that the corresponding additional parameters must be 
input , eg . the floating sun mass and sun center springs . Two 
additional degrees of freedom are added to any system if the 
floating sun parameters are input. 

Phasing constants for mesh time location simulation: 

For a phasing constant of 0 the stiffnesses will be at the pitch 
diameter or time zero, while for a phasing constant of 1, the 
stiffnesses will again be at the pitch diameter but at the end of 
the meshing time. Thus, for phasing constants between 0 and 1, 
the teeth will be located at a corresponding percentage of the 
total mesh time. 


- 51 - 



Other program In puts aff ecting the output: 


The normalized eigenvectors or mode shapes can be printed if 
desired. These will indicate displacements along the lines of 
action for the rotational modes, which are all modes except when 
there is a floating sun. The additional equations for the 
floating sun are for translational movement of the sun center in 
the x and y directions and will lead to translational modes. 

If the full output is requested, the user may find rigid body 
modes that do not always correspond to a natural frequency of 
exactly zero , but are several orders of magnitude less than the 
actual frequencies and are the first frequencies output. This is 
caused by the numerical eigenvalue/eigenvector solution. 

Systems 6,7 and 6 may yield rigid body modes if the stiffnesses 
for the carrier and/or ring gear rim or pin stiffnesses are 
insufficient. For these situations either the carrier and/or the 
ring are acting as rigid bodies. For these systems all the 
frequencies are output, but if the first two frequencies are 
orders of magnitude less than the other frequencies, the pin 
stiffnesses and the results should be carefully examined. 

Input geometry can affect the tooth pair stiffnesses, as well as 
the input torque. 

Finite Element Helical Gear Tooth Analysis: 

This option is initiated by a trigger in location 111. Note that 
the default for helical gear tooth analysis is 0., or the original 
option of 10 independent axial segments, and that no provision is 
made for double helical tooth forms in the finite element 
solution. 

Another point of interest is the numerical solution stability. 
The size of the mesh time step has an effect on the numerical 
solution. This is currently set to be equal to the total mesh 
time divided by 100 for spur gears and the finite element helical 
option. For the original helical option the total mesh time is 
divided by 10. For very low rpms this time step can become 
relatively large and could lead to an unstable solution. 


Flexible Carrier or Flexible Ring Gear Rim Option: 

The flexible carrier/ring gear rim option is triggered by the 
input system type, for either a planetary or star system. The 
additional inputs required are stiffnesses (locations 800-803) In 
addition, the boundary conditions must now be input for each 
individual component (locations 481 to 621). The following 
procedure is recommended to reduce the iterations for convergence: 
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1 . Run the case without a flexible carrier or ring gear rim 
(k ■ 1 or 3). 

2. Estimate the carrier segment or ring gear rim displacement 
(assuming equal segment displacements) from : 

XCO(I) - total output torq ue 

(carrier base radius * pin stiffness) 


or : 


XRO(I) = total output torque 

(ring gear base radius * pin stiffness) 


3 . Using the converged displacement boundary conditions from 

Step 1 , solve for the remaining displacements recalling the 
following tooth pair displacement (along the line of action) 
relationships: 


XOSP(I) - XSO - XPO(I) - XCO(I) 

XORP(I) - XPO(I) - XCO(I) - XRO(I) 

Recall that for planetary systems XRO(I) - 0. and for star 

systems, XCO(I) - 0. 


4. For the initial velocity conditions, one value must be assumed 
(e.g. carrier velocity « 0.) as there are two equations and three 
unknowns . The same relations hold as in Step 3 replacing the 
displacements with velocities. 
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*********************************** 
*********** F178 ****************** 

*********************************** 

* * USER'S MANUAL FOR MULTIPLE MESH* 

** SPUR AND HELICAL GEAR SYSTEMS *** 
**************************** ******* 
********** JUNE 1986 *************** 

*********************************** 

********* PROGRAM RESTRICTIONS ************** 

1. MAXIMUM PLANETS FOR SPUR AND SINGLE HELICAL GEARING 
IS 20 

2. MAXIMUM D LANETS FOR DOUBLE HELICAL GEARING IS 10 

3. MAXIMUM ITERATIONS FOR SOLUTION CONVERGENCE IS 20 

4. MAXIMUM TEETH WITH TOOTH SPACING ERRORS IS S 

5. 10 TEETH ARE CHECKED FOR DYNAMICS DURING TOOTH ERROR PASS 

6. MAXIMUM INVOLUTE CONTACT RATIO IS 3.0 

7. MAXIMUM NUMBER OF SUN TEETH IS 50 FOR RUNOUT OPTION 

*********************************** 

* * *INPUT DATA FOR MULTIPLE MESH SPUR AND HELICAL GEAR PROGRAM * * * 

NOTE: ALL GEOMETRIC INPUT DATA IS IN THE ROTATIONAL PLANE AND 

ALL DEFAULT VALUES ARE 0.0 UNLESS OTHERWISE NOTED. 

ANGLES ARE IN DEGREES 
FORCES ARE IN POUNDS (#) 

LENGTHS ARE IN INCHES (IN) 

MASSES ARE IN (#-SEC**2/IN> 

STRESSES ARE IN (#/IN**2 ) 

TEMPERATURES ARE IN DEGREES FAHRENHEIT 


THE DATA SET REQUIRED TO RUN THE MULTIPLE MESH PROGRAM IS 
DESCRIBED BELOW AND AN EXAMPLE DATA SET WITH CORRESPONDING 
EXAMPLE PROBLEM ARE IN APPENDIX C, 1984 MULTIPLE MESH REPORT NASA 
CR- 174747. A LIST OF GEAR SYSTEM PARAMETERS WITH DESCRIPTIONS, 
LOCATION NUMBERS, AND SOME NONZERO DEFAULT VALUES FOLLOW THE 
INPUT DESCRIPTION. 

CARD Is TITLE CARD— CONTAINS CASE TITLE AND/OR DESCRIPTIVE 
(LINE 1) INFORMATION. AT LEAST ONE COLUMN FROM 

1 TO 40 MUST BE NONBLANK. 

DATA CARDS: CARDS CONTAIN DATA IN THE FORM SHOWN IN FIGURE 1. 

FROM LEFT TO RIGHT THE INPUTS ARE: NUMBER OF DATA 

ITEMS ON CARD, LOCATION NUMBER OF FIRST DATA ITEM 
IN LINE (SEE GEAR SYSTEM PARAMETERS SECTION FOR 
LOCATION NUMBERS AND CORRESPONDING PARAMETER 
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DESCRIPTIONS), AND DATA ITEMS ARE IN 5 FIELDS OF 12 
SPACES. 

CASE TERMINATION CARD: TO TERMINATE A CASE THE LAST LINE MUST 

CONTAIN 0-1. IN COLUMNS 1-4. 

SUBSEQUENT CASES: ANOTHER DATA SET OF THE SAME FORMAT— TITLE 

CARD, DATA CARDS, AND CASE TERMINATION CARD 
—MAY FOLLOW THE CASE TERMINATION CARD. 

JOB TERMINATION CARD: AFTER THE LAST SUBSEQUENT CASE TERMINA- 

TION CARD, A BLANK CARD (LINE) MUST BE 
INCLUDED, OTHERWISE AN ERROR MESSAGE IS 
GENERATED. 

** NOIE ** ALL NUMBERS MUST BE REAL EXCEPT THE NUMBER OF 
ITEMS, THAT IS, A DECIMAL POINT MUST BE INCLUDED. IT IS NOT 
NECESSARY TO INPUT ZERO VALUES AS THEY WILL DEFAULT TO ZERO 
UNLESS OTHERWISE SPECIFIED. 

A BRIEF DESCRIPTION OF OUIPUT FOLLOWS THE GEAR SYSIEM 
PARAMETER SECIION. 

*** THE FOLLOWING EXAMPLE IS FROM THE 1984 NASA CR *** 

IN APPENDIX C THE TITLE CARD READS ’EXAMPLE DATA SET’. THE NEXT 
LINE IN THE DATA SET CONTAINS, FROM LEFT TO RIGHT, S, 
CORRESPONDING TO THE FIVE DATA HEMS ON TO BE ON THAT LINE, THEN 
THE LOCATION NUMBER (1.) CORRESPONDING TO THE FIRST INPUT DATA 
TO BE PUT IN COLUMNS 13-24, FOLLOWED BY THE DATA ITEMS: LEVEL 
(2.), THE DIAMETRAL PITCH (8.4667), THE PRESSURE ANGLE (22.5), 
COAST SIDE PRESSURE ANGLE AND HELIX ANGLE. (BOTH COAST SIDE 
PRESSURE ANGLE AND HELIX ANGLE - 0 FOR STANDARD SPUR GEARS AND 
WOULD DEFAULT TO 0 IF LEFT BLANK) 

MOST OF THE INPUT VARIABLES HAVE SUFFICIENT EXPLANATION IN 
THE GEAR SYSTEM PARAMETERS SECTION, HOWEVER THE FOLLOWING 
PARAMETERS SHOULD BE CALCULATED AS FOLLOWS. 

*** EQUIVALENT MASSES: 

SUN EQUIVALENT MASS - J / (BASE RADIUS OF SUN)**2 

WHERE J IS THE MASS MOMENT OF INERTIA FOR THE SUN. THE OTHER 
COMPONENTS ARE CALCULATED SIMILARLY, USING THE CORRESPONDING 
MOMENTS AND BASE RADII. 

*** PHASING CONSTANTS: 

FOR EQUALLY SPACED PLANETS THE SUN PHASING CONSTANTS ARE 
DETERMINED BY ASSUMING THE FIRST PLANET MESH HAS A PHASING 


I* 
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CONSTANT OF ZERO. THE REMAINING SUN-PLANET PHASING CONSTANTS 
ARE DETERMINED BY: 

(PLANET # - 1 ) *THE FRACTIONAL REMAINDER FROM DIVIDING THE 

NUMBER OF SUN TEETH BY THE TOTAL NUMBER OF PLANETS 
WHERE THE 'FRACTIONAL REMAINDER’ INDICATES THE SPACING DIFFER- 
ENCE BETWEEN PLANETS. 

THE PHASING CONSTANTS FOR THE RING-PLANET MESHES ARE DETERMINED 
THE SAME WAY IF THE PLANET HAS AN ODD NUMBER OF TEETH, USING THE 
NUMBER OF RING TEETH. IF THE PLANETS HAVE AN EVEN NUMBER OF 
TEETH THE CONSTANTS ARE CALCULATED FOR THE SUN MESHES IHEN 0.5 
IS ADDED TO EACH TO OBTAIN THE RING-PLANET MESH PHASING 
CONSTANTS (KRP) , DUE TO THE RING GEAR BEING .5 OFFSET FROM THE 
SUN GEAR. 

EXAMPLE CORRESPONDING TO APPENDIX C: 

SUN PLANET MESH— 14/3 - 12 2/3 

PLANET 1, KSP ( 1 ) - 0. 

PLANET 2. (2-1) *(2/8) - .6666667 
KSP ( 2 ) - .6666667 

PLANET 3, (3-l)*(2/3)-l. 3333333 
KSP(3 ) - .8333333 

RING PLANET MESH— (EVEN NUMBER OF PLANET TEETH) 


PLANET 

1. 

KRP(l) - 

KSP(l) + .5 - .5 

PLANET 

2, 

KSP(2) 
KRP(2) - 

♦ .5 - 1.1666667 
.1666667 

PLANET 

3, 

KSP(3) 

♦ .5 - .8333333 



KRP ( 3 ) - 

.8333333 


*** BOUNDARY CONDITIONS: 

AN INITIAL ESTIMATE FOR DISPLACEMENT BOUNDARY CONDITIONS MAY 
BE OBTAINED BY DIVIDING THE STATIC TOOTH LOAD AT THE PITCH RADIUS 
BY THE TOOTH SPRING RATE AT THE PITCH RADIUS. 
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*Z*7.*7.*Z*Z*Z*Z 1986 UPDATE *7.*Z*Z*Z*Z*7.*Z*Z*Z 

SEE THE 1986 CR REPORT FOR A SUMMARY OF THE LOCATIONS FOR THE 
ADDITIONAL INPUTS REQUIRED FOR THE NEW OPTIONS. THE SAME 
LEVELS APPLY, WHERE THE NEW OPTIONS CAN BE USED WITH ANY LEVEL 
OF INPUT. NOTE THAT THE LOCATIONS 10 TRIGGER THE PLOTS HAVE 
BEEN CHANGED. ALSO NOTE THAT THE NEW OPIIONS THAT ADD DEGREES 
OF FREEDOM ARE SENSITIVE 10 THE INIIIAL BOUNDARY CONDITIONS, 
SEE CR FOR INIIIAL FLEXIBLE CARRIER BOUNDARY CONDITION INPUT 
RECOMMENDATIONS. 


*Z*Z*Z*Z*7.*Z*Z*Z*Z*Z*Z*Z*Z*Z*Z*Z*Z*7.*Z*Z*Z*Z*Z*Z*Z*Z*Z*Z*Z* 


I **** ***** LEVELS OF INPUT 

j LEVEL I INPUT 

THE FIRST INPUT ITEM WILL BE THE LEVEL DESIRED, FOLLOWED BY ITEMS 
2 THROUGH 51 AND WHERE NOTED. THE OTHER VALUES WILL DEFAULT. 

THE MATERIAL PROPERTIES, LOCATIONS 22 TO 27, DEFAULT FOR STEEL. 
THE TOLERANCES ARE SET TO ZERO OR DEFAULT VALUES BELOW, ERRORS 
, ARE SET TO 0. THERE ARE NO PROFILE MODIFICATIONS. OTHER NONZERO 

DEFAULT VALUES ARE: 

1 SURFACE ROUGHNESS - 25 

OIL TEMPERATURE - 180 F 
MATERIAL CONSTANT - 0.0528 
DAMPING - .02 
OIL TYPE MIL-L-23699 

| ANY DESIRED PLOTS CAN BE OBTAINED, SEE LOC 651-658, 

: A CHECK ON INPUT DATA CAN BE MADE, LOC # 699. 

IF CONTACT RATIO IS HIGH (GREATER THAN 2) IT MUST BE INPUT IN LOC 
54 & 55. 

i 

i 

LEVEL II INPUT 

j THE MINIMUM INPUT FOR THIS LEVEL WOULD BE LEVEL 1 DATA PLUS 

DAMPING, FLASH TEMPERATURE DATA, SOLUTION ITERATION DATA, AND 
PHASING CONSTANTS. DEFAULT VALUES ARE ZERO UNLESS INDICATED 
! OTHERWISE BELOW. THE MAXIMUM INPUT WOULD INCLUDE THE ITEMS 

REQUIRED FOR LEVEL 1 PLUS ALL OTHER ITEMS EXCEPT LOCATIONS 64, 

| 65, 91-100, AND 175-192. 

LEVEL III INPUT 

LEVEL 3 REQUIRES ALL ITEMS TO BE INPUT, UNLESS DEFAULT VALUES ARE 
INDICATED BELOW. THE MINIMUM INPUT FOR THIS LEVEL WOULD BE 
LEVEL 2 DATA. 

**** NOTE: FOR HIGH CONTACT RATIO GEARS (CR > 2.), THE USER MUST IN- 

PUT THE CR (LOC 54 AND 55). THE PROGRAM CALCULATES RATIO 
IF CR < 2. 
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*** GEAR SYSTEM PARAMETERS ****************** 


LOC 

NAME 

DESCRIPTION 

* * * 

* * * * 

**************************** 

1 

LEVEL 

TRIGGER FOR LEVEL OF INPUI DATA 

2 

DP 

DIAMETRAL PITCH (NORMAL PLANE) 

3 

PSANG 

SUN-PLANET DRIVE SIDE PRESSURE ANGLE fi PD 
RING-PLANET COAST SIDE PRESSURE ANGLE 8 PD 

i 4 

! 

PSANGB 

SUN-PLANET COAST SIDE PRESSURE ANGLE 9 PD 

RING-PLANET DRIVE SIDE PRESSURE ANGLE 9 PD 

** 0.0 IF PRESSURE ANGLES FOR DRIVE & COAST SIDES EQUAL 

5 

PSIO 

HELIX ANGLE § PD 
**0.0 FOR SPUR GEARS 

6 

N1 

NUMBER OF TEETH ON SUN GEAR 

7 

N2 

PLANET GEARS 

8 

N3 

RING GEAR 

9 

FW1 

AXIAL FACE WIDTH OF SUN GEAR 

10 

FW2 

PLANET GEARS 

11 

FW3 

RING GEAR 

12 

N 

NUMBER OF PLANET GEARS 

13 

K 

IF K-l PLANETARY SYSTEM, I. E. .RING GEAR FIXED 


IF K-2 STAR SYSTEM , I. E., PLANET CARRIER FIXED 
*** IF K-3 DIFFERENTIAL SYSTEM, I. E. NEITHER RING OR CARRIER 
FIXED 

7.7. ** IF K-4 NON PLANETARY, I. E. , NO RING AND NO CARRIER 

SUN-INPUT, PLANET-OUTPUT 

7.7. IF K-S NON PLANETARY, I. E. , NO SUN AND NO CARRIER, 

PLANET INPUT, RING OUTPUT 

* IF K-6 PLANETARY SYSTEM WITH FLEXIBLE PLANET CARRIER 

* IF K-7 STAR SYSTEM WITH FLEXIBLE RING GEAR RIM 

******* K - 6 AND K - 7 CURRENTLY UNSTABLE RESULTS ****** 

* IF K-6 OR K-7 LOC 88 & 89 ARE REQUIRED FOR PLANET 
CARRIER OR RING GEAR RIM STIFFNESS VALUE ALONG THE 
LINE-OF-ACTION AND AN INTERFACE STIFFNESS. 

** FOR K-4 AND RUNOUT ERROR, INPUT LOC 200. 

*** FOR K-3, INPUT LOC 18 -21. 

7.7. FOR K-4 OR 5, THE SINGLE MESH PROGRAM WILL GIVE THE 
RESULTS IN LESS CPU TIME DUE TO AN EXACT SOLUTION. 

ALSO, FOR K - 5, THERE IS A BUG IN THE DEPTH TO MAX 
SHEAR CALCULATION IN THE MULTI -MESH PROGRAM. 

14 TORQ AXIAL INPUT TORQUE ON SUN GEAR (OR PLANET IF K=5) 

15 RPM INITIAL AXIAL ROTATIONAL SPEED OF SUN GEAR 

(OR PLANET IF K=5) 

16 RPMF FINAL AXIAL ROTATIONAL SPEED OF SUN GEAR FOR SPD RANGE 

17 INTVL NUMBER OF MAIN INTERVALS THE SPEED RANGE DIVIDED INTO 

** FOR ONE RPM ONLY, RPM-RPMF AND INTVL-1.** 
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STEP SIZE - ( RPMF - RPM)/INTVL, THIS STEP SIZE WILL 
BE AUTOMATICALLY BY 5 IN AREAS OF PEAK LOADS, THUS 
REFINING THE INCREMENT. 


******* DIFFERENTIAL SYSTEM INPUTS ***************************** 


18 

19 

20 
21 


TOUTC OUTPUT TORQUE FROM CARRIER 

TOUIR OUTPUT TORQUE FROM RING 

RPMC CARRIER RPM 

RPMR RING RPM 


******* GEAR MATERIAL PROPERTIES ******************************* 


22 El 

23 E2 

24 E3 

25 MU1 

26 MU2 

27 MU3 


YOUNGS MODULUS * E-06 OF SUN GEAR (DEFAULT - 30.) 


PLANET GEARS “ “ “ 

RING GEAR “ “ “ 

POISSONS RATIO OF SUN GEAR (DEFAULT - .30) 

PLANET GEARS 

RING GEAR “ “ “ 


******* GEAR EQUIVALENT MASSES ********************************* 

28 MS EQUIVALENT MASS OF SUN GEAR ABOUT ROTATIONAL AXIS 

29 MC PLANET CARRIER 

80 MR RING GEAR 


31-50 MP ( I ) PLANET GEARS 

I - PLANET NUMBER 

*** NOTE: FOR K-6 OR K-7, CARRIER OR RING MASS FOR TOTAL 
UNIT, NOT SEGMENTS. 


******* TRIGGER FOR DOUBLE HELICAL GEARING ********************* 


51 DBHEL IF - 0.0 SINGLE HELICAL GEARING IF PSIO .NE. 0.0 
IF - 1.0 DOUBLE HELICAL GEARING IF PSIO .NE. 0.0 


******************************** 

END OF LEVEL 1 INPUT (UNLESS HIGH CONTACT RATIO) 
******************************** 

******* GEAR MESHING DAMPING RATIOS **************************** 

52 ZSP DAMPING RATIO OF SUN-PLANET MESHES 

53 ZRP “ “ “ RING-PLANET MESHES 

* * * NOTE: IF CONTACT RATIO > 2, CRSP AND CRRP MUST BE INPUT, BUT 

THE PROGRAM WILL CALCULATE FOR CR < 2 ONLY. 

54 CRSP INVOLUTE CONTACT RATIO OF SUN-PLANET MESHES 

55 CRRP “ “ “ “ RING-PLANET MESHES 
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******* ACTIVE FACE WIDTH VALUES ******************************* 


**NOTE: LOC 56 - 59 FOR HELICAL GEARS ONLY, NOT NECESSARY FOR 

FINITE ELEMENT OPTION 


56 

LISP 

INACTIVE 

SUN-PLANET 

FACE 

WIDTH 

ON LEFT 

57 

L2SP 

tt 

a tt 

tt 

(( 

“ RIGHT 

58 

L1RP 

INACTIVE 

RING-PLANET 

FACE 

WIDTH 

ON LEFT 

59 

L2RP 

it 

tt tt 

ft 

tt 

“ RIGHT 


******* SUN-PLANET PROFILE MODIFICATION INPUT DATA ************* 

LOCATIONS 66- 68 AND 75-78 ARE ONLY USED FOR PROFILE MODIFICATION 
TABLES, THEY ARE NOT USED IN THE DYNAMICS. 

60 PCTSOD SD AS A PERCENT OF SOD (Z> 

61 PCTSOE SE AS A PERCENT OF SOE (Z) 

62 DELD DISENGAGEMENT TIP RELIEF (IN.), MINIMUM 

68 DELE ENGAGEMENT TIP RELIEF (IN. ), MINIMUM 
**NOIE: SICN CONVENTION— POSITIVE INPUT, MATERIAL REMOVED 

****** XNVDX AND XNVEX ARE NOT REQUIRED FOR LEVEL 2 * * • * 

64 XNVDX DISENGAGEMENT PROFILE MODIFICATION SHAPE FACTOR 

65 XNVEX ENGAGEMENT 

********************************* 

66 PMODSP PERCENT OF TIP RELIEF ON SUN GEAR, ENGAGEMENT 

67 PMDSPD PERCENT OF IIP RELIEF ON SUN GEAR, DISENGAGEMENT 

- Z IIP RELIEF ON PLANET 

68 DLSP TOLERANCE AT START OF PROFILE MODIFICATION 

78 DLIOL TOTAL IIP RELIEF TOLERANCE (SUN ♦ PLANET 

AND/OR RING + PLANET, DEPENDING ON SYSTEM TYPE) 

SOD - LENGTH OF DISENGAGEMENT PART OF THE LINE OF ACTION 
SOE - “ “ ENGAGEMENT “ “ “ “ “ “ 

SD & SE ARE THE SEGMENTS OF THE SOD & SOE THAT ARE USED, WHEN 
PROFILE MODIFICATIONS ARE MADE. 

******* RING-PLANET PROFILE MODIFICATION INPUT DATA ************ 

69 PTS0D3 SD AS A PERCENT OF SOD 

70 PTS0E3 SE AS A PERCENT OF SOE 

71 DELD3 DISENGAGEMENT TIP RELIEF 

72 DELE3 ENGAGEMENT TIP RELIEF 

***** LOCATION 73 AND 74 NOT REQUIRED FOR LEVEL 2 * * * * * 

73 XNVDX3 DISENGAGEMENT PROFILE MODIFICATION SHAPE FACTOR 

74 XNVEX3 ENGAGEMENT PROFILE MODIFICATION SHAPE FACTOR 
******************************** 

75 PMODRP PERCENT OF TIP RELIEF ON PLANET GEAR, ENGAGEMENT 
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76 PMDRPD PERCENT OF TIP RELIEF ON PLANET GEAR, DISENGAGEMENT 

- % TIP RELIEF ON RING GEAR 

77 DLRP IOLERANCE AT START OF PROFILE MODIFICATION 

78 DLTOL TOTAL TIP RELIEF TOLERANCE, FOR BOTH SUN+PLANEI 

AND/OR RING + PLANET 

****** MESH MODIFICATION DUE TO FACE WIDTH CROWNING *********** 

80 LECSP LENGTH OF FACE WIDTH CROWN OF ENGAGEMENT FOR SP MESH 

81 LDCSP “ “ “ “ “ “ DISENGAGEMENT FOR SP MESH 

82 DLECSP ENGAGEMENT EDGE RELIEF FOR SP MESH 

83 DLDCSP DISENGAGEMENT EDGE RELIEF FOR SP MESH 

84 LECRP LENGTH OF FACE WIDTH CROWN OF ENGAGEMENT FOR RP MESH 

85 LDCRP “ “ “ “ “ “ DISENGAGEMENT FOR RP MESH 

86 DLECRP ENGAGEMENT EDGE RELIEF FOR RP MESH 

87 DLDCRP DISENGAGEMENT EDGE RELIEF FOR RP MESH 


****** TOOTH PAIR COMPLIANCE DATA ***************************** 

90 WIOK IF - 0.0 PREPROCESSOR CALULATES TOOTH PAIR COMPLIANCE 
i IF - 1.0 INPUT COMPLIANCE DATA IN LOCATIONS 91-101 

| (FOR FINITE ELEMENT HELICAL COMPLIANCE, LOC 111) 

***** 91 - 100 NOT REQUIRED FOR LEVEL 2************ 

i 

SINGLE TOOTH PaIR SPRINGRATE OF SUN-PLANET MESHES 
COMPLIANCE CONSTANT <S/S0)**1 OF SUN-PLANET MESHES 

2 

3 

4 

SINGLE TOOTH PAIR SPRINGRATE OF RING-PLANET MESHES 
COMPLIANCE CONSTANT (S/S0)**1 OF RING-PLANET MESHES 

2 

3 

4 

WHERE S - IS THE LENGTH ALONG THE LINE OF ACTION TO THE POINT 
OF ENGAGEMENT 

AND SO - THE LENGTH OF THE LINE OF ACTION 
********************************* 

101 HRTZSP HERTZ STRESS FOR COMPLIANCE CALCULATION FOR SUN-PLANET 

102 HRTZRP RING-PLANET 

** THE PROGRAM WILL USE 101 AND 102 AS CONSTANTS FOR THE COM- 
PLIANCE CALCULATIONS IF INPUT. DEFAULT WILL CALCULATE 


91 SPKSP 

92 COMASP 

93 COMBSP 

94 COMCSP 

95 COMDSP 

96 SPKRP 

97 COMARP 

98 COMBRP 

99 COMCRP 
100 COMDRP 



THE STATIC LOAD/PITCH RADIUS FOR THE CONSTANTS 

103 XPSPSP IF = 0.0 PLANE STRESS IS ASSUMED FOR SUN-PLANET 

IF - 1.0 PLANE STRAIN IS ASSUMED FOR SUN-PLANET 

104 XPSPRP IF - 0.0 STRESS RING-PLANET 

IF - 1.0 STRAIN RING-PLANET 

111 CONVEX IF = 0.0 HELICAL TOOTH IS DIVIDED INTO 10 INDEPENDENT 

AXIAL SEGMENTS 

IF - 1.0 FINITE ELEMENT ROUTINES USED TO GENERATE AN 
EQUIVALENT SPUR GEAR COMPLIANCE CURVE 
*£&****£ FLOATING SUN PARAMETERS ******************************* 


115 

KFSX 

LINEAR 

SPRING 

IN 

THE 

X 

DIRECTION 

AT 

SUN 

CENTER 

116 

KFSY 

LINEAR 

SPRING 

IN 

THE 

Y 

DIRECTION 

AT 

SUN 

CENTER 

117 

DFSX 

LINEAR 

DAMPER 

IN 

THE 

X 

DIRECTION 

AT 

SUN 

CENTER 

118 

DFSY 

LINEAR 

DAMPER 

IN 

THE 

Y 

DIRECTION 

AT 

SUN 

CENTER 

119 

FSMS 

ACTUAL 

SUN GEAR MASS 

(NOT EQUIVALENT MASS) 


******* TOOTH PAIR GEOMETRIC DATA ****************************** 


120 BYPASS IF - 0.0 GEOMETRIC PREPROCESSOR IS USED 

IF - 1.0 GEOMETRIC DATA MUST BE INPUT IN LOCATIONS 
121-132 


121 

RADIOS 

MAX. RADIUS TO BASE OF FILLET OF SUN GEAR 

(ROOT 

122 

RADTOP 

PLANET GEARS 

RADIUS) 

123 

RADTOR 

RING GEAR 


124 

RADTIS 

MAX. RADIUS 10 TIP OF TOOTH OF SUN GEAR 


125 

RADTIP 

PLANET GEARS 


126 

RADTIR 

RING GEAR 


127 

TOOTHS 

NOMINAL TOOTH THICKNESS AT PD OF SUN GEAR 


128 

TOOTHP 

PLANET GEARS 


129 

TOOTHR 

RING GEAR 


130 

RADFIS 

MAX. FILLET RADIUS OF SUN GEAR 


131 

RADFIP 

PLANET GEARS 


132 

RADFIR 

RING GEAR 



FLASH TEMPERATURE DATA ********************************* 

140 

SROU 

SURFACE ROUGHNESS (RMS, MICROINCHES) 


141 

TOIL 

OIL INLET TEMPERATURE 


142 

CFMAT 

MATERIAL CONSTANT 


143 

OILTYP 

TYPE OF OIL IN GEARBOX 



1 FOR MIL-L-23G99 

2 FOR MIL-L- 7S0S 

3 TO BE DETERMINED 
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4 TO BE DETERMINED 
(0 DEFAULTS TO TYPE 1) 


******* SOLUTION ITERATION DATA ******************************** 



150 TOLER ITERATION TOLERANCE FOR BOUNDARY CONDITION CONVERGENCE 

(Z * 100, DEFAULT FOR LEVEL 1 - .01) 

151 NLOOP NUMBER OF ITERATIONS FOR BOUNDARY CONDITION CONVERGENCE 

NLOOP MAXIMUM - 20, DEFAULT - 20. 

****** TOOTH PAIR GEOMETRIC TOLERANCE DATA ******************** 

160 DLROS TOLERANCE ON TIP RADIUS OF SUN GEAR (DEFAULT - .002) 

161 DLROP PLANET GEARS 

162 DLROR RING GEAR “ 

163 BSE EDGE BREAK ON IOPLAND OF ALL GEARS (DEFAULT = .010) 

164 DLRR ROOT RADIUS TOLERANCE OF ALL GEARS (DEFAULT - .005) 

165 DLCDN CENTER DISTANCE TOLERANCE (TOWARD TIGHT MESH) ALL GEARS 

166 DLCDP CENTER DISTANCE TOLERANCE (AWAY FROM TIGHT MESH) 

167 DLMBT MACHINE BACKLASH TOLERANCE ALL GEARS (DEFAULT « .002) 

NOTE: THE TOOTH ROOT FILLET RADIUS TOLERANCE - .005 IN THE PREPROCESSOR 
******* 8-DIMENSIONAL FACTORS FOR SPUR GEARS ONLY ************** 
***** 175 - 192 NOT REQUIRED FOR LEVEL 2************ 

175 SPRNBS IF - 0. RIM BENDING EFFECTS NOT INCLUDED IN SUN GEAR 

(DEFAULT - 10.E*10) 

IF > 0. BENDING SPRINGRATE ( IN-#/RAD ) OF SUN RIM 

176 SPRNBP IF - 0. RIM BENDING EFFECTS NOT INCLUDED IN PLANET GEAR 

(DEFAULT - 10.E+10) 

IF > 0. BENDING SPRINGRATE ( IN-#/RAD ) OF PLANET RIM 

177 SPRNBR IF - 0. RIM BENDING EFFECTS NOT INCLUDED IN RING GEAR 

(DEFAULT - 10.E+10) 

IF > 0. BENDING SPRINGRATE! IN-#/RAD) OF RING RIM 

178 SPRNLS RADIUS FROM RIM TO PD (REQUIRED WITH SPRNBS) SUN 

179 SPRNLP RADIUS FROM RIM TO PD (REQUIRED WITH SPRNBP) PLANET 

180 SPRNLR RADIUS FROM RIM TO PD (REQUIRED WITH SPRNBR) PLANET 

181 EFWD1S EFFECTIVE SUN FACE WIDTH FACTOR AT TOOTH TIP 

IF - 0. (DEFAULT - 1.0) 

182 EFWD2S EFFECTIVE SUN FACE WIDTH FACTOR AT FILLET 

IF - 0. (DEFAULT - 1.0) 

183 EFWD1P EFFECTIVE PLANET FACE WIDTH FACTOR AT TOOTH TIP 

IF - 0. (DEFAULT - 1.0) 

184 EFWD2P EFFECTIVE PLANET FACE WIDTH FACTOR AT FILLET 

IF « 0. (DEFAULT ** 1.0) 

185 EFWD1R EFFECTIVE RING FACE WIDTH FACTOR AT TOOTH TIP 

IF - 0. (DEFAULT - 1.0) 


- 64 - 


I 


I 



186 

187 

188 

189 

190 

191 

192 


EFWD2R EFFECTIVE RING FACE WIDTH FACTOR AT FILLET 

IF ■= 0. (DEFAULT = 1.0) 

STRSES STRESS DISTRIB. FACTOR FOR END EFFTS AT SUN WIDTH EDGE 
IF - 0. (DEFAULT - 1.0) 

STRSCS STRESS DISTRIB. FACTOR FOR END EFFTS AT SUN WIDTH CENT 
IF - 0. (DEFAULT - 1.0) 

STRSEP STRESS DISTRIB. FACIOR FOR END EFFTS AT PLAN WDTH EDGE 
IF - 0. (DEFAULT - 1.0) 

STRSCP STRESS DISTRIB. FACTOR FOR END EFFTS AT PLAN WDTH CENI 
IF - 0. (DEFAULT - 1.0) 

SIRSER STRESS DISTRIB. FACTOR FOR END EFFTS AT RING WDTH EDGE 
IF - 0. (DEFAULT - 1.0) 

STRSCR STRESS DISTRIB. FACTOR FOR END EFFIS AT RING WDTH CENT 
IF - 0. (DEFAULT - 1.0) 


******* HELIX ANGLE ERRORS ************************************* 


195 

DELPSP 

HELIX ANI 

196 

DELPRP 

« 1 

197 

DELPS2 

IF 

DBHEL 



IF 

DBHEL 

198 

DELPR2 

IF 

DBHEL 



IF 

DBHEL 


“ “ RING-PLANET MESH 

0.0 DELPS2 - 0.0 

1.0 DELPS2 - HELIX ANGLE ERROR OF RIGHT 
HALF OF SUN-PLANET DOUBLE HELICAL GEARS 
0.0 DELPR2 - 0.0 

1.0 DELPR2 - HELIX ANGLE ERROR OF RIGHT 
HALF OF RING-PLANET DOUBLE HELICAL GEARS 


******* TOOTH PAIR SPACING ERRORS ****************************** 


200 DR SUN RUNOUT ERROR FOR EXTERNAL-EXTERNAL SINGLE MESH ONLY 
(K-4), DISPLACEMENT ERROR OF SUN CENTER 

** NOTE: A WRITE (7 ,5029) STATEMENT NEEDS A COMMENT STATEMENT 

REMOVED FOR TIME AND LOADS TO BE OUTPUT, OTHERWISE 
CARDS ARE PUNCHED WHEN A TAPE IS WRITTEN. IT IS IN 
SUBROUTINE STRESS. 


**NOTE: FOR DOUBLE HELICAL GEARS (DBHEL - 1.0) THE FIRST 10 

VALUES OF I IN THE E(I,J) ARRAYS ARE FOR THE LEFT HALF 
OF THE DOUBLE HELICAL GEARS AND THE LAST 10 LOCATIONS 
ARE FOR THE RIGHT HALF OFTHE DOUBLE HELICAL GEARS, 

WHERE I - THE PLANET NUMBER. 

THE ERRORS SHOULD BE PUT ON THE SECOND OR THIRD TOOTH TO ILLUSTRATE 
THE DYNAMICS. 


221 TO 240 ESP (1,1 )SUN-PLANET ERROR ARRAY FOR TOOTH 1 

241 TO 260 ESP (1,2) 2 
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2G1 TO 2S0 ESP (I, 3) 3 
281 TO 300 ESP(I,4> 4 
301 TO 320 ESP ( 1,5) 5 

321 TO 340 ERP(I,1)RING-PLANET ERROR ARRAY FOR TOOTH 1 
341 10 360 ERP(I,2) 2 
361 TO 380 ERP(I,3) 3 
381 TO 400 ERP(I,4) 4 
401 TO 420 ERP(I,5) 5 


******* PLANET GEARS PHASING CONSTANTS ************************* 

421 TO 440 KSP(I) SUN-PLANET PHASING CONSTANTS ARRAY 
441 TO 460 KRP(I) RING-PLANET PHASING CONSTANTS ARRAY 

******* INITIAL BOUNDARY CONDITIONS **************************** 

481 TO 500 XOSP(I) SUN-PLANET INITIAL DISPLACEMENT BOUNDARY COND 

OR SUN DISPLACEMENT FOR K-6 OR K-7. 

501 TO 520 XORP(I) RING-PLANET INITIAL DISPLACEMENT BOUNDARY COND 

OR PLANET DISPLACEMENT FOR K-6 OR K-7. 

521 TO 540 X1SP(I) SUN-PLANET INITIAL VELOCITY BOUNDARY COND 

OR SUN VELOCITY FOR K-6 OR K-7. 

541 TO 560 XIRP(I) RING-PLANET INITIAL VELOCITY BOUNDARY COND 

OR PLANET VELOCITY FOR K-6 OR K-7. 

FLOATING SUN INITIAL CONDITIONS 

561 CARRIER DISPLACEMENT 

581 CARRIER VELOCITY 

601 RING DISPLACEMENT 

621 RING VELOCITY 

641 SUN CENTER DISPLACEMENT IN X DIRECTION 

642 SUN CENTER DISPLACEMENT IN Y DIRECTION 

643 SUN CENTER VELOCITY IN X DIRECTION 

644 SUN CENTER VELOCITY IN Y DIRECTION 

******* PROGRAM PLOT SELECTIONS FOR SPUR GEARS ONLY ************ 

651 PLTLD IF - 0.0 

IF - 1.0 LOAD PLOTS, NORMALIZED LOAD — DYNAMIC/STATIC 

652 PLTPV IF - 0.0 

IF - 1.0 PV PLOTS, PRESSURE SLIDING VELOCITY 

653 PLTHS IF - 0.0 

IF - 1.0 HERTZ STRESS PLOTS 

654 PLTFT IF - 0.0 

IF = 1.0 FLASH TEMPERATURE PLOTS 

655 PLTSS IF = 0.0 

IF - 1.0 SUN GEAR HEYWOOD STRESS PLOTS 

656 PLTPSS IF - 0.0 
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IF - 1.0 SUN-PLANET PLANET GEAR HEYW00D STRESS PLOTS 
657 PLTPRS IF - 0.0 

IF - 1.0 RING-PLANET PLANET GEAR HEYWOOD STRESS PLOTS 
65S PLTRPS IF - 0.0 

IF - 1.0 RING GEAR HEYWOOD STRESS PLOTS 

******* PROGRAM CHECK RUN ************************************** 

699 CHECK IF - 0.0 REGULAR RUN TO COMPLETION 

IF - 1.0 EXITS PROGRAM BEFORE DYNAMIC SOLUTION 

TO ALLOW CHECKING OF INPUT DATA AND THE 
GEOMETRIC PREPROCESSOR RESULTS 


*********************************** 
****** FLEXIBLE PLANET CARRIER OR RING GEAR RIM DATA ********** 

800 CARRK1 IF K-6, PLANET CARRIER STIFFNESS ALONG-THE-LINE-OF- 

ACTION (#/IN) 

801 CARRK2 IF K-7 RING GEAR RIM STIFFNESS ALONG-THE-LINE-OF- 

ACTION (#/IN) 

802 KFLX1 PLANET/CARRIER ’PIN' OR INTERFACE STIFFNESS 

803 KFLX2 PLANET/RING ’PIN’ OR INTERFACE STIFFNESS 

********************* NATURAL FREQUENCIES OPTION ***************** 

805 NFREQ TRIGGER FOR NATURAL FREQUENCY OPTION 

- 1, OUTPUTS FREQUENCIES IN TERMS OF INPUT RPM 

- 2, OUTPUTS FREQUENCIES IN TERMS OF INPUT RPM 

AS WELL AS THE FREQUENCIES WITH EIGENVECTORS 

** SEE USER’S MANUAL IN 1986 CR FOR ADDITIONAL SYSTEM 
TYPES THAT ARE AVAILABLE 
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***** OUTPUT DESCRIPTION 


THE FOLLOWING LIST INDICATES THE INFORMATION AND RESULTS THAT 
MAY APPEAR IN THE OUTPUT, IN THE ORDER THEY WILL APPEAR. 

CASE TITLE— -THE TITLE AND/OR OTHER INFORMATION INPUT ON THE FIRST 
DATA CARD. 

TIP MODIFICATION — IF THERE IS INSUFFICIENT TIP CLEARANCE, AN 
ADJUSTMENT IS MADE INTERNALLY AND A MESSAGE OUTPUT UNTIL 
SUFFICIENT CLEARANCE IF OBTAINED. 

INVOLUTE MODIFICATION TABLE— SHOWS THE PROCESSED RESULTS OF 
ANY INPUT PROFILE MODIFICATIONS. TABLES ARE PRINTED FOR 
ENGAGEMENT AND DISENGAGEMENT WHICH INCLUDE MINIMUM AND 
MAXIMUM INVOLUTE MODIFICATIONS, DIAMETER AND CORRESPONDING 
ROLL ANGLE FOR THE MODIFIED PORTION OF THE PROFILE. THIS 
WILL BE OUTPUT FOR THE SUN-PLANET MESH AND/OR RING-PLANET 
MESH IN CONJUNCTION WITH THE CORRESPONDING INPUT DATA. 

INPUT DATA— THE INPUT DATA AND PREPROCESSED GEOMETRIC DATA IS 
PRINTED FOR SUN-PLANET MESH AND/OR RING-PLANET MESH. 

FLEXIBILITY— IF THE FLEXIBLE PLANET CARRIER OPTION IS IN EFFECT 
A MESSAGE APPEARS THAI INDICATES THIS. 

ADDITIONAL INPUT DATA— NUMBER OF PLANETS, EQUIVALENT MASSES, ETC. 

COMPLIANCE CONSTANTS— CALCULATED CONSTANTS FOR FOURTH ORDER 
COMPLIANCE EQUATION. 

ITERATED BOUNDARY CONDITIONS— BOUNDARY CONDITIONS PRINTED, FOL- 
LOWED BY THE CURRENT SPEED BEING EXAMINED AND CORRESPONDING 
MAXIMUM LOADS (FOR SUN-PLANET AND/OR RING-PLANET MESHES). 

IF A SPEED SURVEY WAS RUN, THE SPEED CORRESPONDING TO THE 
OVERALL MAXIMUM LOAD FROM THE RANGE CALCULATED IS OUTPUT 
WITH THE MAXIMUM LOAD. THIS IS THE SPEED USED FOR THE 
REMAINING CALCULATIONS, I.E. STRESS. 

MAXIMUM VALUES— TABLE(S) OF MAXIMUM VALUES CALCULATED ARE PRINTED 
FOR EACH SUN-PLANET MESH AND/OR RING-PLANET MESH. THESE 
WILL APPEAR FOR THE NO ERROR SOLUTION AND EACH ERROR SOLU- 
TION. TWO ABBREVIATIONS APPEAR — PV-PRESSURE SLIDING VELO- 
CITY AND PD-PITCH DIAMETER. 

ERROR MATRIX— IF TOOTH SPACING ERRORS ARE INPUT, OR GENERATED 
FOR RUNOUT SOLUTION, A TABLE OF ERRORS IS PRINTED BEFORE 
THE TABLES OF MAXIMUM VALUES. 

PLOTS — WILL BE FOR EACH PLANET MESH. 
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IF THERE ARE TOOTH PAIR SPACING ERRORS INCLUDED .PLOTS ARE 
GENERATED FOR THE NO ERROR CASE ONLY , I.E. THERE ARE NO 
PLOTS GENERATED FOR THE ERROR CASES. 
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APPENDIX B: FIRST ORDER DIFFERENTIAL EQUATIONS 


A. Floating Sun Equations: 


The equations can be transformed to first order equations as 
follows. 


Let; : 



-1 C I d 

m 


N 


s i=l ■»! ““ “1 T ± i x ^s Pl sln “i 


sin « A + i L sin 
i-1 ‘ 

_d * x i " k x x a ] 


(2a) 

(2b) 


[ Ji d ®Pi 008 °i 


N 

I L cos a. 
i-1 s *i 1 


+ d 7 X 3 


k y x 4 3 


(3a) 

(3b) 


N 

t I d 
m 


N 

+ l 


c ±ri « P± *s Pi l£l r P± ^r P± 
N 

I X* + _ 

i-1 ? Pi 


N N 

+ L s Pi + .1 + ] 


- X c 


(4a) 

(4b) 


N 

t I d. 


N 

+ I K 


i=i r Pi ir Pi iti r P 


+ F. 3 


(5a) 

(5b) 
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B. Flexible Carrier Equations 

The carrier and ring gear equations are reduced to first order 
equations as follows. 


Let 


'10 


■'ll 


■'ll 


(6a) 

(6b) 

(6c) 


x °ir sj V ( V x Pii' ~ £ W\i~ \r 


+ K c (X C - X n ) +K 


'C sPi " < Lrp 1 ~°^o ^ ( " x - - + x * - > 


"i-i 


- '« w \ 


V 


Let 


(7) 

(carrier) 


X - y 

10 Yr l 

(8a) 

X ” V 

11 Yr l 

(8b) 

X - Y 

r ll r i 

(8c) 


r il m r dr Pl CX Pn X °n X » 1 x ) a Lr Pi 


Z 1 k r. , ( - x rj . + X P . )US, (X, 

r i 


n *i_i T. ' “r. ' ^ a r. '• A r " x r > 
r i-l 1 0 1 0 *r, 1 % i+1 o 


out r • -r. - “b r 


X„ - X „ 

10 11 


(9) 

(ring) 
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The equation for the torque constraint reduces as follows. 
Let 

(10a) 
(10b) 
(10c) 
(10d) 

* ~ (T °ut 7 \ - j^oTo (11) 



APPENDIX C: FORTRAN LISTING 


The program listing is available from NASA Project Manager or 
Contractor upon request. 
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APPENDIX D : NOMENCLATURE 


A,B,C,D 

C 


d 


rPi 

l sp ± 


d. 


' rp ij 

' Sp ij 


f I 

I 

[K] 

K c 

c i 

K 


K. 


k 


rp i 


k 


sp. 


k 


k 


gear tooth coefficients 

single tooth pair compliance (in/ lb) 

single tooth pair compliance at pitch radius 
(Reciprocal of spring rate) (in/lb) 

ring-planet tooth pair damping (lb s/in) 

sun-planet tooth pair damping (lb s/in) 

damping at sun center In x-direction (lb s/in) 

damping at sun center in y-direction (lb s/in) 

tooth spacing error, ring-planet mesh 

tooth spacing error, sun-planet mesh 

gear mesh frequencies 

integer multiplier of gear mesh frequencies 

stiffness matrix (lb/in) 

planet /carrier pin stiffness (lb/in) 

carrier segment stiffness (lb/ in) 

planet/ring stiffness (lb/in) 

ring gear rim segment stiffness (lb/in) 

ring-planet tooth pair stiffness (lb/ in) 

sun-planet tooth pair stiffness (lb/in) 

spring rate at sun center in x-direction (lb/in) 
spring rate at sun center in y-direction (lb/in) 
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rp± 

L sp ± 

[H] 

H 

m 

m c 

m r 

N 


r 

rpm 

S 


X LOA d 

XN 

y 

Y rp. 


£ sp ± 


- ring-planet tooth pair load for planet mesh i 

- sun-planet tooth pair load for planet mesh i 

= mass matrix, a diagonal matrix (lb.) 

- artificial mass for total torque constraint 

*» translational mass of the sun gear (lbs s**2/in) 

- rotational (equivalent) mass of the planet 

carrier (lbs s**2/in) 

■ rotational (equivalent) mass of the ring gear 
(lbs s**2/in) 

- number of planets 

- Carrier Base Radius (in) 

- Ring Gear Base Radius (in) 

- speed of driving gear 

- motion along line of action from pitch line 

" Reference distance along line of contact 
for tooth pair dompliance coefficients 

* displacement of the sun center in x-direction 
(in) 

- Sun-planet displacment along the line of action 

due to the floating sun (in) 

- Number of teeth on driver 

- displacement of the sun center in v-direction 

(in) 

= ring-planet tooth pair displacement for planet mesh i 

= sun-planet tooth pair displacement for planet 
mesh i 

= 0 O + '•'i ~ $ (planetary) 

= - § - 0 r (star) 
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1 

1 

rp ji 

- 

Cam modification, ring-planet mesh 

0 __. 

Sp ji 

- 

Cam modification, sun-planet mesh 

; ^ 

- 

Ring-Planet tooth pair spring rates 

^SPi 

- 

Sun-Planet tooth pair spring rates 

e c 

- 

carrier angle of rotation 

e r 

- 

ring angle of rotation 

X 

= 

2 

O) 

T out 

- 

Output Torque 

$ 

- 

pressure angle 

$ 

. rp ji 

- 

Identify function for ring-planet tooth pair 
contact 

! V 

- 

Identify function for sun-planet tooth pair 
contact 

! *i 

- 

relative angular position of planets 
2 tt (i-l)/N , i - 1,2,3, . . .N 

! U) 

1 

- 

natural frequency (rad/s) 

! x 

sp ji 

1 

Subscri 

pts 

0 or 1 depending on whether the tooth contact is 
on the profile modification cam or not 

1 

i 

: 0 


carrier 

! 

1 

= 

planet mesh 

J 

= 

tooth pair mesh 

P 

1 

= 

planet gear 

I 

- 

ring gear 

s 

= 

sun gear 

i 
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